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ABSTRACT 

We report a periodicity of ~1 day in the highly elevated X-ray emission 
from the protostar V1647 Ori during its two recent multiple- year outbursts of 
mass accretion. This periodicity is indicative of protostcllar rotation at near- 
breakup speed. Modeling of the phased X-ray light curve indicates the high- 
temperature (~50 MK), X-ray-emitting plasma, which is most likely heated 
by accretion- induced magnetic reconnection, resides in dense (>5xl0 10 cm -3 ), 
pancake-shaped magnetic footprints where the accretion stream feeds the new- 
born star. The sustained X-ray periodicity of V1647 Ori demonstrates that such 
protostellar magnetospheric accretion configurations can be stable over timescales 
of years. 
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1. Introduction 

Strong X-ray emission and collimated jets from newborn stars, so-called protostars, 
indicate that energetic magnetic activity plays an important role in star formation. However, 
the thick gaseous envelopes that obscure protostars at visual wavelengths complicate efforts 
to probe their innermost active regions. Low-mass stars like the Sun in the formation phase 
(hereafter, young stellar objects or YSOs) gradually accumulate mass from the parent cloud 
before igniting nuclear fusion in the stellar core. Matter falling from the cloud first forms 
a disk around the central star and matter in the innermost disk gradually accretes onto it. 
The young star’s gravitational potential can accelerate infalling matter up to a few hundred 
km s -1 . When the accreting material collides with the stellar surface, it is shock-heated to 
a few MK; this thermalized plasma can emit soft (E < 1 keV) X-rays (Kastner et al. 2002; 
Sacco et al. 2010). 

Matter cannot fall from the disk onto the central star unless it loses a significant fraction 
of its angular momentum. Coupling of the large scale magnetic fields of the central star to 
those of the innermost disk is suspected to prompt the momentum transfer (e.g., Shu et al. 
1996; Hayashi et al. 1996). Magnetic reconnections triggered by the interactions may be 
involved in the ejection of a fraction of the infalling disk matter, along with most of the 
angular momentum, out of the system; the rest is accreted onto the central star. Bipolar jets 
or collimated outflows seen ubiquitously outside of the YSO’s envelopes are likely launched 
by magnetohydrodynamic or magnetocentrifugal processes (e.g., Reipurth & Bally 2001, and 
references therein). 

Although protostellar magnetospheric accretion models are widely accepted, and have 
gained support from observations of specific, relatively evolved objects (e.g., classical T Tauri 
stars, Donati et al. 2011a, b), the accretion geometry - and hence the validity and applicabil- 
ity of such models at earlier protostellar evolutionary stages - remains uncertain in the case 
of younger, more deeply embedded objects (Johns-Krull et al. 2009). X-ray observations of 
rapidly accreting objects may offer a means to probe this magnetospheric accretion process. 
The bulk of a star’s mass is accreted during very early (observationally, Class 0 and I) phases, 
when the star is still deeply embedded in its parent cloud. Such protostars tend to emit hard 
(>2 keV) X-rays with occasional rapid (<1 day) flares (Imanishi et al. 2001). Some of the 
flares have been proposed as arising in large magnetic loops that connect the inner accre- 
tion disk and the stellar surface (i.e. , in a non-solar- type magnetic structure), which may be 
loaded by accreting material (Tsuboi et al. 1998; Montmerle et al. 2000; Favata et al. 2005). 
However, systematic studies of X-ray-emitting pre-main sequence stars in the Orion nebula 
suggest that the occurrence of the largest magnetic loops, associated with major flare events, 
is not dependent on the presence of circumstellar disks (Getrnan et al. 2008a, b; Aarnio et al. 


3 


2010 ). 

Over the last decade, a small number of YSOs have been found to display sudden in- 
creases in mass accretion rate, by as much as a few orders of magnitude (Hartmann & Kenyon 
1996). Such eruptive YSOs are crudely classified as either FU Ori (FUor) or EX Lupi (EXor) 
types, wherein the former generally display major outbursts over timescales of decades, and 
the latter generally display smaller, shorter-duration, outbursts. One such eruptive YSO, 
V1647 Ori — which does not fit neatly into either the FUor or EXor eruptive class (see discus- 
sion in Teets et al. 2011) — went into outburst during the period 2004-2006 and then, again, 
in 2008 (the latter eruption is evidently still ongoing). This Class I YSO (Muzerolle et al. 
2005, Principe et al. 2012, in prep.) offered the first direct evidence that X-ray activity 
surges with increases in mass accretion activity (Kastner et al. 2004, 2006; Teets et al. 2011). 
During these eruptions, the X-ray luminosity of V1647 Ori increased by two orders of magni- 
tude and the plasma temperature reached ~50 MK (Kastner et al. 2004; Grosso et al. 2005; 
Grosso 2006; Kastner et al. 2006; Hamaguchi et al. 2010; Teets et al. 2011). Subsequently, 
other YSOs that experienced mass accretion outbursts were found to have similarly enhanced 
or strong X-ray activity (Audard et al. 2010; Skinner et al. 2006, 2009; Grosso et al. 2010). 
The gravitational potential of YSOs is not deep enough to thermalize plasma to such a 
high temperature; yet the elevated levels of X-ray emission observed from YSOs undergoing 
accretion outbursts (such as V1647 Ori) strongly indicate an intimate connection between 
accretion activity and X-ray emission in these objects. In order to determine where in the 
star/disk system the hard X-rays actually originate, we must ultimately understand the 
origin of the high energy activity during accretion outbursts. To this end, we have reana- 
lyzed all X-ray observations of V1647 Ori in outburst, in search of temporal evidence that 
might reveal the site(s) and, perhaps, mechanism(s) responsible for its enhanced high-energy 
emission during these events. 


2. Data Sets 

Since the onset of the first outburst in 2004, V1647 Ori has been monitored 17 times 
with three major X-ray observatories: 13 times with Chandra (Weisskopf et al. 2002), three 
times with XMM-Newton (Jansen et al. 2001) and once with Suzaku, (Mitsuda et al. 2007). 
Six Chandra observations during the first outburst and an XMM-Newton observation during 
the second outburst did not collect enough photons for timing analysis. Table 1 lists the 
9 data sets used for the timing analysis and the 2 additional data sets added to the folded 
light curve plots (Figs. 1, 2). Hereafter, individual Chandra , XMM-Newton and Suzaku 
observations are designated CXO, XMM and SUZ, respectively, subscripted with the year, 
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month and day of the observation. 

We reprocessed data with the following calibration versions: DS 7.6 or later for the 
Chandra data, SAS 10.0.0 or later for the XMM-Newton data and version 2.2.11.22 for 
the Suzaku data. The basic analysis follows procedures described in earlier papers - 
Chandra (Kastner et al. 2006; Teets et al. 2011), XMM-Newton (Grosso et al. 2005), and 
Suzaku (Hamaguchi et al. 2010). A notable departure from those procedures is that we 
converted event arrival times to the solar barycenteric time system, which differs by up to 
300 sec from the terrestrial time. We generated background subtracted light curves for pho- 
tons with energies in the range of 1—8 keV, binned into 2 ksec (= AT) intervals, using the 
HEAsoft 1 analysis package. We developed Python codes for the cross-correlation, chi-square 
and physical model studies. 


3. Period Search 

We identified strong similarities in eleven separate X-ray light curve observations of 
V1647 Ori obtained during the two outbursts with Chandra , XMM-Newton and Suzaku. 
These flux variations — an order of magnitude on timescales of hours — superficially resemble 
stellar magnetic flares, but there is reason to doubt this interpretation, given the spectral vari- 
ation and frequency of flux rises (Grosso et al. 2005; Kastner et al. 2006; Hamaguchi et al. 
2010). Fig. 1 displays all 11 light curves with time offset and flux normalization according 
to the x 2 study described below. We first focus on the light curve with the longest duration, 
obtained with XMM-Newton in 2005 (Fig. 1, top ; the numbering scheme in the following sen- 
tence corresponds with the labels at the top of the figure). The X-ray flux i ) stays constant 
for ~20 ksec, ii ) rises by a factor of ~5 in ~15 ksec, Hi) keeps an elevated level for 30 ksec 
with marginal spikes and dips, and iv ) falls gradually to the original flux level on a similar 
time scale. The XMM-Newton light curve in 2004 apparently matches with i) — iii), and the 
Suzaku light curve in 2008 starts from iv) and connects to i)—iii). Although the Chandra 
light curves (Fig. 1 bottom) are of more limited durations and have lower photon statistics, 
each of those light curves also matches with one or more parts of profile i) — iv). 

A cross-correlation analysis provides quantitative support to these qualitative similari- 
ties (see Appendix A for details). The XMM-Newton light curves in 2004 and 2005 correlate 
strongly (0.92) when they align at their observation starts and the former is shifted back- 
ward by 6 ksec. The XMM-Newton light curve in 2005 and the Suzaku light curve in 2008 
correlate strongly (0.82) when they are folded by a period of 86 ksec and the Suzaku light 


1 http: / /heasarc. gsfc.nasa.gov /docs/ software /lheasoft / 
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curve shifts backward by 30 ksec. 

Given the similarities of these light curves, we search for a more accurate period and set 
of phase shifts that match both the shapes and the timings of all the available light curves (see 
Appendix B for details). First, we generate a template light curve with a 1.23 day span by 
combining the XMM-Newton and Suzaku light curves that are shifted in time and normalized 
in flux based on the cross-correlation study. We repeat the template light curve with a 
frequency below 1.45 day” 1 , giving a phase offset between 0.0— 1.0 and a flux normalization 
for each light curve to account for the long-term variation (Teets et al. 2011). Based on our 
preliminary analysis, we also introduce a phase gap between the first and second outbursts. 
We fit the template to all light curves with the least y 2 method and derive the minimum y 2 
value of 317.9 [y 2 /d.o.f = 1.88 ( d.o.f . = 169)] at /o = 0.98929 day” 1 (P =87.3 ksec) and 
&ga P = —0.383. The best-fit period is close to the rotation period of the Earth; however, 
because none of the (space-based) observatories whose data are analyzed here obtain data 
on a daily cadence, this period cannot be an artifact of our observing protocol. We therefore 
conclude that V1647 Ori displays periodic variation of its X-ray emission with a period of 
~1 day. 

One obvious potential origin for the X-rays would be the rotation into and out of our 
line of sight of a localized region of X-ray plasma on the central star. Rises and falls in 
the light curves would correspond to appearances and disappearances of the localized X-ray 
bright spot. The flux transitions take ~20% of one cycle (Fig. 3, Appendix C), suggesting 
that the spot has a significant size or height compared with the radius of the central star. 

We assume a uniform circular spot with a tip-cut cone shape and simulate an X-ray light 
curve for each combination of the spot radius, height, latitude and stellar inclination angle. 
We find that no single spot can produce both the low flux interval (</> ~0. 0—0.2, 0.8— 1.0) and 
the high flux interval (0 ~0.4— 0.6), so we adopt a bipolar geometry and add a fainter spot 
at the opposite longitude and latitude of the brighter spot whose shape is identical to that 
of the bright spot (see Appendix C for details). The best-fit result (y 2 =98.2, d.o.f. =36) is 
obtained with a bright spot to faint spot brightness ratio of 5, with the spots having radii 
of 0.32 A* and heights of 0.01 R *, and with the bright spot found at a stellar latitude of 
~ —49°. The stellar inclination — the tilt of the polar axis toward our line of sight — is 
~68°. This model, shown in Fig. 3, reproduces the low and high flux levels and the rise and 
fall of the represented light curve well. The possible dip at = 0.5 may be reproduced by a 
partial eclipse of the bright spot by its own accretion flow or the disappearance of the faint 
spot behind the central star if it has a smaller size and higher latitude than those assumed 
in the fit above. Excesses at $ = 0.4 and 0.8 after the flux rise and fall may represent 
asymmetries of the spot shapes or the presence of additional smaller spots. The best-fit for 
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the inclination angle of the rotation axis is similar to the angle (~61°) estimated from an 
infrared light echo study (Acosta-Pulido et al. 2007). 

The one-cycle light curve can also be reproduced given plasma configurations that differ 
somewhat from the geometry described in Fig. 3. For example, the faint hot spot — which is 
required in the preceding model so as to reproduce the low flux levels during phases intervals 
$ = 0.0 — 0.2, 0.8 — 1.0 — may shift toward the latitude and longitude directions with 
respect to the bright hot spot, or can be replaced by a constantly visible plasma in the form, 
e.g., of a halo around the star. Furthermore, the bright hot spot could really be a complex 
of multiple, smaller hot spots, instead of a uniform single spot. However, in any plasma 
configuration, the flux increase during the phase interval $ = 0.2 — 0.8 constrains the bright 
hot spot (or the envelope of multiple hot spots) to have the approximate size, height and 
latitude described above. 


4. Discussion 

The confidence ranges of the stellar inclination and the latitude of the bright spot (Fig. 4 
left) exclude solutions involving a bright spot in the hemisphere facing us, a bright spot at 
a high latitude, and a pole-on view of the system. This result is consistent with modeling 
of strong fluorescent iron lines observed in Suzaku and Chandra spectra (Hamaguchi et al. 
2010; Teets et al. 2011), which suggests that a significant fraction of hard X-ray-emitting 
plasma is occulted from view. The spot radius can be as large as ~0.5i?*, while its height 
should be lower than ~0.1i?* (Fig. 4 right). The spot — shaped like a thin, extensive plate 
- is similar to shape to the mass accretion footpoints of neutron stars, white dwarfs, and 
the Earth’s aurorae. 

No periodic variation such as that seen in the X-ray regime has been identified in 
optical or infrared observations of V1647 Ori, most likely because the optical and infrared 
emission mostly comes from the disk. Based on a pre-outburst bolometric luminosity and 
stellar effective temperature, Aspin et al. (2008) roughly estimated the mass and radius of 
V1647 Ori at ~O.8M 0 and ~5 R Q , respectively. The stellar radius is slightly larger than the 
distance at which an orbiting body around a O.8M 0 star would have an orbital period of 
about one day. 2 Thus, the central star must be rotating at a speed close to break-up velocity 
(i.e. , rotating at the Keplerian velocity at the stellar radius). 


2 Hereafter, we assume a stellar radius of 4 solar radii — the approximate maximum radius that a 0.8 Mq 
star with rotation period ~1 day can maintain without breaking up. 
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Fig. 4 right also plots the product of the plasma density squared ( n 2 ), the plasma 
filling factor ( 77 ), and the cube of the stellar radius (i?*), using the plasma emission measure 
determined during the Suzaku observation in 2008. There is no solution at log ijn 2 (i2* / 4i ? 0 ) 3 
<2x21 cm -6 . Since 77 ^1 and i?* < 4 R@, we estimate n >5xl0 10 cnr 3 . Being only a 
lower limit, this result for density — which is similar to those of the densest active stellar 
coronae (e.g., Ness et al. 2004) — may place V1647 Ori in the density regime inferred for 
the footpoints of free-fall accretion (as measured for isolated T Tauri stars via line ratios of 
helium-like ions; e.g., Kastner et al. 2002 for TW Hya; see also Porquet et al. 2010). The 
magnetic held B should be stronger than B >100 Gauss at the footpoints, in order to confine 
such a dense hot plasma (i.e., if the magnetic pressure is stronger than the plasma pressure; 
plasma /3 = nkT / (B 2 / 8tt) < 1). 

The phase gap that we infer between the first and second outbursts may be caused 
by drift of the magnetic poles on the stellar surface or a change in the stellar rotational 
frequency. For the latter case, we can replace the phase gap in our model ephemeris with a 
frequency derivative. In doing so, we find a similarly good solution at a similar frequency 
with a small derivative (see Appendix B). 

The observed X-ray variation of V1647 Ori can be naturally explained by rotational 
modulation of X-ray bright spots. Coronal active regions can produce such rotational X- 
ray modulation (Flaccomio et al. 2005). However, earlier studies (Kastner et al. 2004, 2006; 
Teets et al. 2011) indicate that large increases in X-ray flux from V1647 Ori during the out- 
bursts are very closely correlated with the (accretion-driven) optical/infrared flux; therefore, 
these X-ray eruptions are most likely accretion- driven as well. The geometrical model de- 
scribed in §3 (Fig. 3) supports such a model, in that it indicates that the X-ray-emitting 
plasma lies at or very near the foot-points of mass accretion streams at the stellar surface. 

Since the profile of the X-ray light curves did not change remarkably between the ob- 
servations, the intrinsic X-ray luminosity of the hot spots, that is, the accretion-induced 
magnetic reconnection activity, varies on a timescale of a week or longer. Given the large 
overall variation in the amplitude of the X-ray flux from outburst to outburst (and even 
during outbursts) (see Fig. 2), it is evident that the large-scale magnetic held configuration 
of the V1647 Ori star/disk system is preserved, even as the protostar undergoes dramatic 
changes in accretion rate. We suggest two possible mechanisms that may generate this con- 
dition: (*) the mass accretion how stably disrupts the magnetic helds at the footpoint (e.g., 
Brickhouse et al. 2010), or (ii) the rotational shear between the star and the disk continu- 
ously twists the stellar bipolar magnetic helds (e.g., Goodson et al. 1997, see also Fig. 5). 

More evolved YSOs also have faint, hard X-ray emission from hot plasma. This emission 
has usually been explained as due to a blend of emission from multiple micro-hares (e.g., 


Caramazza et al. 2007). Our result demonstrates that the mass accretion activity also can 
generate hot plasma constantly by sustained magnetic reconnection. The same mechanism 
may also operate on those YSOs with weaker mass accretion activity. 


5. Conclusion 

We have discovered rotational modulation of X-ray emission from the Class I protostar 
VI 647 Ori via analysis of 11 X-ray light curves obtained with the Chandra , XMM-Newton 
and Suzaku observatories during this YSO’s two recent mass accretion outbursts. Based on a 
cross-correlation study and period search, we determined a rotational period of ~1 day, with 
either a phase gap apparent between the two outbursts or frequency variation. The single- 
cycle light curve can be successfully reproduced by emission from two hot spots on opposite 
poles on the stellar surface. The star rotates with a period of ~1 day, close to the break-up 
velocity for a 0.8 M 0 star with a radius of ~4 R & . The hot spots likely cover significant 
fractions of the stellar surface, while the height (~0.01 72* ) may be negligible compared with 
the stellar radius. The hot spot size and the plasma emission measure indicate relatively 
high plasma density (>5xl0 10 cm -3 ), also pointing to an origin in accretion hot spots. This 
result clearly demonstrates that hard X-ray-emitting plasma can be present in long-lived 
accretion footprints at the surfaces of protostars, and thereby constrains the geometry of 
magnetospheric accretion in early (Class I) protostellar evolutionary stages. 


This work is performed while K.H. was supported by the NASA’s Astrobiology Institute 
(RTOP 344-53-51) to the Goddard Center for Astrobiology (Michael J. Mumma, P. I.). J.K.’s 
research on X-rays from erupting YSOs is supported by NASA/GSFC XMM-Newton Guest 
Observer grant NNX09AC11G to RIT. This research has made use of data obtained from 
the High Energy Astrophysics Science Archive Research Center (HEASARC), provided by 
NASA’s Goddard Space Flight Center. 


Facilities: Chandra (ACIS), XMM-Newton (EPIC), Suzaku (XIS) 


A. Cross-Correlations of the Light Curves 


We cross-correlate light curves of XMM 040404 , XMM 0503 24 , and SUZ 08100 8, which have 
good photon statistics. We define the cross-correlation index r, 


where Xi and y t are count rates of the i - th bins of two light curves, d is a delay in units of 
2 ksec time bins, x and y are averages of bins that contributes to the cross-correlation. The 
index r ranges between —1 and 1. Two light curves are identical when r = 1. 

We first cross-correlate XMM 0 40404 {xf) with XMM050324 {yf) (Fig. 6 left). There is a 
strong correlation of r = 0.92 at d = 3 (6 ksec). 

Since SUZ 0810 o8 and XMM 050 324 have similar observing durations, the number of bins 
that contributes to the cross-correlation becomes smaller at larger delays. The r index 
fluctuates strongly in these regions. We thus require the number of contributing bins to 
be at least 10. When SUZ 08 ioo8 (xf) is shifted backward against XMM050324 {yi), there is 
a strong correlation of r = 0.91 at d = 28 (56 ksec). When XMM050324 {%%) is shifted 
backward against SUZ 08 ioo8 (Vi), there is a strong correlation of r = 0.81 at d — 15 (30 ksec). 
This result indicates that these two light curves folded at a certain period also match well. 
We thus fold both light curves by 40—46 bins (80—92 ksec), average overlapped bins with 
weighted mean values, and cross-correlate them. We found that the r index is a maximum 
of 0.82 when these light curves are folded by 43 bins (86 ksec) and XMM 050 324 {xf) is shifted 
backward by d = 30 (60 ksec) against SUZ 081008 {yf) (Fig. 6 right). 


The cross-correlation tests the similarity of two light curves, but it does not consider 
the time interval between them. We therefore search for an ephemeris that satisfies both 
shapes and timings of all the X-ray light curves, including those detected with Chandra. 

We first generate a template light curve from the XMM 040 404 , XMM 050 324 and SUZ 08 ioo8 
light curves. Based on the cross-correlation study in Appendix A, we shift XMM 040404 and 
XMM 050324 backward by 18 and 15 bins (36 and 30 ksec) against SUZosioos, respectively. 
We normalize these light curves based on the average count rates of overlapped bins and 
average them at weighted mean values. After artificially increasing the time resolution of 
the averaged light curve by 100 times with linear interpolation, we smooth it with a Gaussian 


J2[(xi ~ x)(yi- d - y)} 


r = 



(Al) 


B. Search for the Best Ephemeris with the y 2 Test 
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function with a = 2 ksec so as to minimize statistical noise. The resulting template light 
curve [L ternp (ti)] spans a time period of 1.23 day. 

We define the ephemeris as: 

$(T) = <h 0 + f 0 (T - T 0 ) + ® gap H(T - T gap ) (Bl) 

where /o is frequency, T 0 is the time origin, <h 0 is the phase at T = T 0 , Q gap is the phase gap, 
H is the unit step function and T gap is the time of the phase gap. We then assign phases, 
'k(fj), to the template light curve according to 


'h(ti) = f 0 U 


(B2) 


and discard bins at v h(t i ) > 1. 

We sum bins of the template light curve within the corresponding phase range of each 
bin of a measured light curve, that is, 


I temp ( P ) ^ ^ L bemp ( tj ) 


j=m 


where m and n satisfy 


&(Ti) = $(P) - L$(P)J 
Ad> 

T(t m ) A $'(Ti) - — < f(l ra+1 ) 
*(<*) S *'(Ti) + — < >t((„ + i) 


(B3) 


(B4) 

(B5) 

(B6) 


where |_^J is the floor function and A$ = foAT. We then normalize ItempiTi) for each 
measured light curve by the normalization factor n obs that gives a minimum y 2 value, and 
derive a y 2 value for each set of $ 0 , & gap , and f 0 . 


T^obs — 


X = 


L obs {Tj) I te mp(Ti) / AL obs (Tj)' 

EhempiT^/AL^iT^ 


EE 


obs 


-LobsiTi) Tl obs 1 temp (2"i) 

AL obs (Ti) 


(B7) 

(B8) 


We search for the minimum y 2 value in the ranges of 0 A <f> 0 , & ga p < 1 and /o A 0.82 day -1 . 
We find a minimum y 2 of 317.93 [y 2 /d.o.f. = 1.88 ( d.o.f . =169)] at / 0 = 0.98929 day -1 
(P =87.3 ksec) and & gap = —0.383. The best ephemeris is hence expressed as 


$(T) = 0.98929 (T - 13453.09075) - 0.383 H(T - 14300) 


(B9) 
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where T is the truncated Julian date. In this formula, we re-define the $ origin at an 
intermediate point between the fall and rise timings of the single cycle light curve (see 
Appendix C). 

The phase gap is empirical. The phase change may instead be explained by a frequency 
variation. We hence modify equation (Bl) to allow for a frequency derivative. 

<h(T) = <h 0 + / 0 (T-T 0 ) + / 1 (T-T 0 ) 2 / 2 (BIO) 

We find a similarly good fit of y 2 = 317.22 at / 0 =0.98932 day -1 and /i = — 1.89xl0” 6 day” 2 
at To = 14747.6758097 day. Figs. 7 and 8 show light curves folded with this ephemeris. 

We generate a single cycle light curve from XMM 0 40404 , XMM 050324 and SUZosioos- We 
fold these light curves with equation (B9) and bin all light curves with AT =1984.9 second, 
so that one whole light curve requires exactly 44 bins. We normalize and average these light 
curves in a manner identical to that which produced the template light curve. 


C. Physical Model 


We consider an X-ray point source sitting at longitude 9 , latitude f and height h from 
the stellar surface, viewed from the rotation axis at inclination angle i (see Fig. 9 left). We 
define the 9 origin as the longitude that crosses the opposite side of the central star from 
the Earth at $ = 0. The point source appears in view at a longitudinal difference A 9 e from 
this (0=0) reference point (Figs. 9 right ) that satisfies: 


cosA9 e 


cost sinf + 


' _R*_' 

' Rt+h > 


sini cosf 


The normalized X-ray light curve of the point source is: 


T point ( < I>) 


0 (0 A S < <3> e ) 

1 (4> e AS Al-$ e ) 
0 (1 - <h e < S < 1) 


(Cl) 


(C2) 

(C3) 

(C4) 


where S = $ + A-L$ + Aj anc i <|> e = 27 tA 9 e . 

X-ray spectra of V1647 Ori indicate that the X-ray plasma is optically thin (Kastner et al. 
2006; Grosso et al. 2005; Hamaguchi et al. 2010; Kastner et al. 2004; Teets et al. 2011). This 
means that any portion of the X-ray plasma with a finite size can be seen once the portion 
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emerges from the stellar rim, such that the observed flux at phase $ is the integrated emission 
from the plasma emerging above the stellar rim, i.e., 





(<f), 9, (£>, h, i)S(9, (f>, h)(R* + h) 2 cos(j)dOd(j)dh 


(C5) 


where S(9,(f>,h ) is the X-ray source distribution of the plasma. 


For simplicity, we assume a uniform plasma with a conical shape with angular radius 
r* and height h\ standing upside down on the stellar surface, such that the base position 
of the cone is ’’exposed” (Fig. 10). It is observed at unit flux when in view, such that the 
source distribution is described as 


S(9,4>,h ) = 1/V ( inside ) 
= 0 ( outside ) 


(C6) 

(C7) 


where V is the plasma volume. We define the longitude and the latitude of the cone axis as 
6' and (ft 1 , respectively. A narrow latitudinal strip at ft + Aft (|A0'| < r*) has a half width 
A 9 S that satisfies 


cosA9 s 


cos r* — sin{ft + A ft)sinft 
cos(ft + A ft)cosft 


(C8) 


We numerically integrate equation (C5) for this plasma, i.e., 

r h ' f*' +r * r d ' +A(>a F p oint($, 9 , ft h, i)(R , + hfcosj ) 

/ 0 Jfi-r* Jd'-Ads 


A' C one(^ > j 9 , (j) , r* , h , l) — 


1/ 


dJ9d(j>iE 9) 


The low and high flux phases in the single cycle light curve cannot be reproduced by 
any single spot. We therefore assume two spots with identical shapes in a dipole geometry, 
sitting on opposite sides of the central star. The X-ray light curve is, then, expressed as 

F = g b F cone ($,9',ft,r*,ti,i) + g f F cone ($,9' + Tr,-ft,r*,ti,i) (CIO) 

where g b and gj are un-occulted fluxes of the bright and faint spots, respectively. The y 2 
value is a minimum (98.2) when ft = —49°, r* = 18°, h! =0.01i?*, i =68° and gf/g b = 0.20. 
In this fit, we fixed the longitude 9' at 0, considering the definition of the $ origin. Fig. 3 
illustrates this best-fit model. 


The plasma emitting volume is, 

V e = TjV 


— q I (if* + h) 2 Qdh 

J 0 

= r *) h/ [3R l + 3RM + h ^ 


(Cll) 

(C12) 


3 


(C13) 


where g and fl are the plasma filling factor and the solid angle of the cone, respectively. To 
constrain the plasma density, we combine the standard relation for plasma emission measure, 
E.M. = n 2 V e ( n : plasma density), with equation (C13), 



27r(4i? 0 ) 3 (l - cos r*)h*(l + K + hi/ 3) 


E.M. 


(C14) 


where h * = h' /R*. We calculate the right side of the equation for each combination of r* 
and h* and E.M. = 1.9 x 10 54 cm~ 3 during the Suzaku observation in 2008 (Hamaguchi et ah 
2010). Fig. 4 right shows contours of values for this parameter. 
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Table 1. Analyzed Data Sets 


Abbreviation 

Observation ID 

Start Date 

Exposure 

(ksec) 

Duration 

(ksec) 

Net Count 
(counts) 

First Outburst: 

CXO040307 

5307 

2004 Mar. 7 

5.5 

5.6 

60 

CXOo40322 1 

5308 

2004 Mar. 22 

4.9 

5.0 

10 

XMM040404 

0164560201 

2004 Apr. 3 

29.1 

37.0 

1321 

XMMo50324 

0301600101 

2005 Mar. 24 

79.2 

89.7 

1557 

CXO050411 

5382 

2005 Apr. 11 

18.2 

18.4 

85 

Second Outburst: 

CXO 080918 

9915 

2008 Sept. 18 

19.9 

20.2 

455 

SUZ 081008 

903005010 

2008 Oct. 8 

40.4 

77.5 

1275 

CXO 081127 

10763, 8585 

2008 Nov. 27 

20.0, 28.5 

77.5 

197,143 

CXO 090123 

9916 

2009 Jan. 23 

18.4 

18.6 

240 

CXO 090421 

9917 

2009 Apr. 21 

29.8 

30.2 

258 

XMMl00228t 

0601960201 

2010 Feb. 28 

33.5 

34.0 

163 


Note. — Abbreviation: CXO — Chandra , XMM — XMM-Newton , SUZ — - Suzaku. f The 
datasets are not used for the timing analysis. Net count: Background subtracted photon 
counts between 1—8 keV. Photon counts of all the available instruments are summed for 
XMM-Newton and Suzaku. 
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Fig. 1. — Light curves of V1647 Ori between 1—8 keV obtained with the XMM-Newton , 
Suzaku (top) and Chandra (bottom) observatories. These light curves are folded and normal- 
ized according to the best fit ephemeris (/o =0.98929 day -1 , $ gap = —0.383; see Appendix B). 
The upper and lower horizontal axes in each panel show the rotational phase and time in 
second from the phase origin, respectively. The phase origin is set at the middle of the flux 
fall and rise. Light curves in thin colors are repeats of thick ones in the same colors. Points 
with dotted error bars are not used for the best ephemeris search. Solid purple and green 
lines in the top and bottom panels are the template light curves for the x 2 analysis. Curves 
with narrow widths indicate phase intervals that have been disregarded, as per the algorithm 
described in Appendix B. The arrows at the top depicts the phases of variation mentioned 
in Section 3. 



Observed Flux (10 13 erg s 
10 150 5 


-18- 


u III 


iv 


u m 


IV 


0 


1 1 1 r 

- XMM-Newton 
Suzaku 


0.5 

— i — 


Phase <X> 
1 


-i 1 1 r 


T 


“i 1 1 r 


1.5 

— i — 


— i 1 1 1 — 

— 2004 Apr. 3 ( XMM) 

— 2005 Mar. 24 ( XMM) 

— 2008 Oct. 8 ( Suzaku ) 
2010 Feb. 28 ( XMM) 


lO 



+ 






t h 


++ + -H4+ 

L 


rW: 


++ ■ i 






L 




hU + . 


i A 



ji 


++ i i 1 i 



Fig. 2. — Same as Fig. 1, but with axes in energy flux units. 
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Fig. 3. — Single cycle light curve created from the XMM-Newton and Suzaku light curves 
overlaid with the best-fit result of the symmetric bipolar spot model (red: total, green: bright 
spot, blue: faint spot). The pictures on the top depict locations of hot spots at corresponding 
phases. See Appendix C for the numbers and letters in the label. The subscripts b and f 
stand for the bright and faint spots, respectively. 
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Fig. 4. — Confidence ranges of parameter values obtained from model fitting (Fig. 3) - 
left: the stellar viewing inclination vs. the latitude of the bright spot; right: the radius 
vs. height of the hot spots, in stellar radii. In both panels, the black, red and green 
contours show confidence levels of 66% (Ay 2 =2.3), 90% (Ay 2 =4.6) and 99% (Ay 2 =9.2), 
respectively. The blue stars show the best- fit result. The right panel also shows contours 
of log gn 2 (i£* / 4i? 0 ) 3 = 21, 22, 23, and 24 in solid orange lines as obtained from modeling 
Suzaku data (Hamaguchi et ah 2010, see Appendix C for the derivation.) 
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Fig. 5. — A possible geometry of the V1647 Ori system (mechanism ii). Differential rotation 
between the star and the disk shears the stellar bipolar magnetic fields, and the magnetic 
fields twist and continuously reconnect. Matter, accelerated to v ~2000 km s _1 by the 
magnetic reconnection, collides with the stellar surface, thermalizes to kT ~4 keV (T ~50 
MK ) and emits hard X-rays (See e.g., the equation (9) of Calvet & Gullbring 1998). The 
opposite magnetic pole may have lower mass accretion and emit weaker X-rays. 
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Fig. 6. — Cross-correlation of XMM050324 with XMM040404 {left) and SUZ 0 8ioos {right) when 
these light curves are folded by 43 bins (86 ksec). 
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Fig. 7. — Normalized light curves (top: XMM-Newton , Suzaku, bottom : Chandra ) between 
1—8 keV folded with the ephemeris in the equation (BIO), which assumes a frequency deriva- 
tive instead of a phase gap. See the caption of Fig. 1 for the other details. 
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Fig. 8. — Same as Fig. 7, but with axes in energy flux units. 
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Fig. 9. — Definitions of the coordinate system, 6, <f, h and i (left) and of A 6 e when h = 0 
( right top : edge-on view, right bottom : pole-on view). 
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Fig. 10. — Definition of the parameters 6', ft, h! , Aft, Ah and A 6 S for the conically shaped 
plasma model. 


